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ABSTRACT 

We have discovered 21 Rotating Radio Transients (RRATs) in data from the Green Bank Telescope 
(GBT) 350-MHz Drift-scan and the Green Bank North Celestial Cap pulsar surveys using a new can¬ 
didate sifting algorithm. RRATs are pulsars with sporadic emission that are detected through their 
bright single pulses rather than Fourier domain searches. We have developed RRATtrap, a single-pulse 
sifting algorithm that can be integrated into pulsar survey data analysis pipelines in order to find 
RRATs and Fast Radio Bursts. We have conducted follow-up observations of our newly discovered 
sources at several radio frequencies using the GBT and Low Frequency Array (LOFAR), yielding 
improved positions and measurements of their periods, dispersion measures, and burst rates, as well 
as phase-coherent timing solutions for four of them. The new RRATs have dispersion measures 
(DMs) ranging from 15 to 97 pccm“^, periods of 240 ms to 3.4 s, and estimated burst rates of 20 to 
400 pulses hr“® at 350 MHz. We use this new sample of RRATs to perform statistical comparisons 
between RRATs and canonical pulsars in order to shed light on the relationship between the two pop¬ 
ulations. We find that the DM and spatial distributions of the RRATs agree with those of the pulsars 
found in the same survey. We find evidence that slower pulsars (i.e. P > 200 ms) are preferentially 
more likely to emit bright single pulses than are faster pulsars (P < 200 ms), although this conclusion 
is tentative. Our results are consistent with the proposed link between RRATs, transient pulsars, and 
canonical pulsars as sources in various parts of the pulse activity spectrum. 

Subject headings: pulsars: general - surveys - methods: data analysis 
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1. INTRODUCTION 

Rotating Radio Transients (RRATs) are a r ecently dis¬ 
covered class of radio-emitting neutron stars (McLaugh- 
lin et al.|2006), which are characterized by sporadic emis¬ 
sion ot individually detectable single pulses. Unlike reg¬ 
ular pulsars, RRATs typically emit a bright pulse, then 
have no detectable emission for many rotations, before 
another such pulse is emitted. The known RRAT popula¬ 
tion is still very small, with only ^100 RRATs knowip^ 
of which only ^25 have phase-connected timing solu¬ 
tions, making the characterization of this source class 
difficult. 

The reason for RRATs’ irregular emission is not yet 
understood, although various models ex ist. _ One pro¬ 

posed explanation by Weltevrede et al.| (2006) is that 
RRATs are simply faint pulsars with high pulse-to-pulse 
amplitude variability, such that only the few brightest 
pulses from these sources are detected. Another sugges¬ 
tion is that RRAT emission is extrinsically modulated, 
fo r example by an asteroid be lt ([Cordes fc Shannon|2008 
Li] 2006). Zhang et al. (2007) alternatively propose that 
RRATs may be old neutron stars whose emission is start¬ 
ing to turn off, or that they are related to nu lling pulsars. 
Specifically, Burke-Spolaor & Bailes (2010) suggest that 
RRATs are an extreme case of nulling pulsars. It is also 
possible that RRATs are explained by some combination 
of the above scenarios. 

Nulling pulsars are intermittent sources that typically 
emit several consecutive pulses, followed by a short pe¬ 
riod of no pulses, until the next set of pulses is detected. 
In this case, most RRATs would be pulsars with nulling 
fractions of > 99%. It may also be that the nulli ng frac¬ 


tion of pnUars increases ove r their lifetime (e.g. |Ritch- 
ings||1976t |Wang et al.||2007|) . RRATs would thus repre 
sent a late evolutionary phase in a pulsar’s life according 
to this hypothesis. This relation between RRATs and 
nulling pulsars is supported by the discovery of a source 
that switches modes between a bright pulsar with a low 
nulling fract ion and a classic RRAT (Burke-Spolaor & 
Bailesj[2010|). 


However, the link between RRATs, canonical pulsars, 
and other apparent pulsar subclasses is not yet clear. It 
has been proposed that RRATs, some of which have ex¬ 
hibited long periods and magnetic fields somewhat higher 
than the average slow pulsar field, may be related to mag- 
netars and X-ray dim isolated neutron star s (XDINS) 


(e-g- 

McLaughlin et al. 

2006 

et al. 



Popov 


class of pulsars, or just at the tail of the distribution 
in terms of pulse-to-pulse variability, the enormous bias 
that exists against finding them in pulsar surveys argues 
that they may represent a large fraction of the radio- 
detectable neutron star (NS) popula tion, perhaps even 
outnu mbering regular radio pulsars (iMcLanghlin et al 


2006). Moreover, Keane & Kramer ( |2008| find that 


when the radio pulsar population is corrected to include 
the potentially large RRAT contribution, the Galactic 
core-collapse supernova rate is too low to explain the NS 
birth rate. This problem is solved by those models that 
state that RRATs are evolutionarily linked to other pul- 

See the ‘RRATalog’ for details on all known RRATs: http: 
//astro.phys.wvu.edu/rratalog 


sar classes. If, on the other hand, RRATs do represent a 
distinct class of sources, this problem remains unsolved. 
Finding and studying more RRATs is therefore an im¬ 
portant step in understanding their nature, and better 
understanding the pulsar population as a whole. 

Along with RRATs, over the last few years there have 
been discoveries of other fast (sub-second timescales) 
transients, characterized by isolated bright dispersed sin- 
gle pulses. The first of these was the “Lorimer burst” 
(Lorimer et al. 2007), a pulse with a seemingly extra- 
galactic origin due to its high dispersion measure (DM), 
though its extragalac tic origin has been questioned (e.g. 
Kulkarni et al. 2014). Another phenomenon is “pery- 
tons,” which are short frequency-swept pulses at high 
DMs, but are thoug ht to be associated with te r restrial 
atmospheric events (Burke-Spolaor et al. 2011). |Keane 


et al. ([2010|) reported another high-UlVI "burst, though 
its Dlvl excess was not sufficiently large to determine 
whether it was cosmological or a burst from a distant 
Galactic RRAT. Finally, in the last two years, a phe¬ 
nomenon known as Fast Radio Bursts (FRBs) emerged, 
which are bright, high-DM pulses that are plausibly ex¬ 
tragalactic and associated with a catacly smic event based 
on their non-repeating n ature (see, e.g., Thornton et al. 


20131 |Spitler et aD[2014|). Moreover, all of the fast tran- 
ab 


sients described above have only been found at a fre¬ 
quency of ~ 1.4 GHz so far. Gharacterizing new tran¬ 
sient, highly dispersed millisecond radio bursts found in 
single-pulse searches is therefore important. 

Here we report on a search for RRATs in two radio pul¬ 
sar surveys conducted at 350 MHz, the Green Bank Tele¬ 
scope Drift-scan survey (hereafter ‘Drift-scan’) and the 
Green Bank North Gelestial Cap (GBNCC) survey, using 
a new search technique. This resulted in the discovery of 
21 new RRATs. We have followed up these RRATs us¬ 
ing the Green Bank Telescope (GBT) and Low Frequency 
Array (LOFAR) telescopes, resulting in improved posi¬ 
tions, periods and DMs. We have also obtained phase- 
coherent timing solutions for four of the RRATs. 

We describe the two surveys in Section[2 , the data anal¬ 
ysis, RRATtrap candidate sifting algorithn ^ and discov¬ 
eries in Section!^ and the follow-up analysis in Section 
We discuss our results in the context of the general pulsar 
population in Section and conclude in Section 

2. SURVEY OBSERVATIONS 

The Drift-scan and GBNCC surveys were both con¬ 
ducted using the Robert C. Byrd Green Bank Telescope, 
whose 100-m diameter, unblocked aperture, and low- 
noise receivers make it an extremely sensitive and pow¬ 
erful tool for radio astronomy. Here we describe these 
surveys, including their sky coverage, sensitivities, and 
pulsar discoveries. We summarize observational param¬ 
eters for the two surveys in T able ^ More detai ls on 
the Drift-scan are presented in [Boyles et a l.| (|2Q13|) and 
Lynch et al.j (|2013|), and on the GBNCC inlbtovall et al. 
(|2014]. - - - 

2.1. Green Bank Teleseope Drift-sean survey 

The GBT Drift-scan survey was a radio pulsar survey 
conducted during the northern summer of 2007, while 

RRATtrap is available at http://github.com/ckarako/ 
rrattrap 
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the azimuth track of the telescope was undergoing re¬ 
pairs. Although the GBT is normally fully steerable, 
this survey took advantage of the time otherwise spent 
idle by pointing the telescope at several fixed elevations, 
and collecting data as the sky drifted over it. Overall, 
this survey produced 134 TB of data over 1491 observ¬ 
ing hours, and covered about 10300 square degrees of the 
sky. The coverage of the Drift-scan spanned declination 
(6) ranges of -8° < (5 < +38° and -21° < 6 < +38°, 
depending on the telescope’s azimuth (which varied be¬ 
tween two values), and is shown in Figure 

Although the recording of data was continuous for each 
observing session, the data were divided into pseudo¬ 
pointings, each 140 s in duration, roughly corresponding 
to the time it takes a point on the sky to pass through 
the full-width half-maximum (FWHM) of the telescope 
beam. This resulted in about 30000 pseudo-pointings, 
each independently analyzed as described in Section 
Each pointing was searched up to a DM of 1015 pc cm , 
with ^10000 DM trials. The RRATtrap algorithm was 
run on the single-pulse results from all beams, constitut¬ 
ing 1491 hr X 0.28 deg^ = 418 hr deg^ of sky. Note that 
the Drift-scan analysis was conducted before the FRB 
phenomeno n was unambiguously confirmed (Thornton 
et al.|2013|), and thus the data were not se arch ed to such 
high DlVls as for the GBNCC (see Section |2.2|. A repro¬ 
cessing of the data with a search of higher DMs is now 
beginning. 

This survey yielded 35 new pulsars. Amongst the 
new pulsars are some exotic systems, including a mil¬ 
lisecond pulsar in a hierarchical triple system with two 
white dwar f companions (PSR J0337+1715; Ransom 
et al.||2014|), and a radio pulsar/low mass X-ray binary 
trans ition object (PSR J1023-1-0038; Archibald et al. 
2009). Furthermore, applying the new single-pulse sift- 


mg techniques developed here, we have discovered 11 
RRATs, and 5 as-yet-unconfirmed RRAT candidates. 
These RRAT discoveries are pre sente d in Table and 
elaborated upon in Sections |3.2| and |4.4| The RRAT 
candidates, as well as the number and duration of re¬ 
observations, are presented in Table 

2.2. Green Bank North Celestial Cap survey 


The GBNCC survey is an ongoing pulsar survey that 
began in 2009 and which also uses the GBT at 350 MHz. 
Unlike the Drift-scan, this survey was designed with the 
goal of uniformly covering the sky visible to the GBT. 
The survey makes use of the Gree n Bank Ultimate Pul sar 
Processing Instrument (GUPPI; |DuPlain et al.||2008|), a 
powerful backend that is now routinely used for pulsar 
observations at the GBT. 

The overall region to be observed is divided into point¬ 
ings, and over the search campaign, each pointing is ob¬ 
served for 120 s. The region that will be covered at the 
completion of the survey is the entire GBT-visible sky, 
—45° < S < 90°, and consists of about 125000 point¬ 
ings. Figure shows the current sky coverage of the 
survey. So far, over 66000 pointings have been observed, 
representing about half of the total planned survey area. 
In the early processing, each pointing was searched up 
to a DM of 500 pccm“^, with ^17500 DM trials, but 
in the last year the maximum DM has been increased 
to 3000 pccm“^ in order to search for highly-dispersed, 
possibly extragalactic bursts. The previously analyzed 


pointings have not yet been reprocessed with this in¬ 
creased maximum DM. This DM threshold was chosen 
to be substantially larger than the highest obse rved DM 
ofFRBs, ^1100 pccm“^ ( Thornton et al.|2013 ), and was 
limited to 3000 pccm“^, since at this DM the time delay 
across the GBNCC observing band is ^60 s, which is half 
of the observation duration. Approximately 200 hr deg^ 
have so far been analyzed with RRATtrap. 

It is estimated that the GBNCC survey will be about 
2.5 times more sensitive to low-DM pulsars at high Galac¬ 
tic latitudes compared to pas t and some ongoing pulsar 
surveys of the same region (Stovall et al. 2014). The 
data analysis for this survey consists of the tecnniques 
described in Section along with improved methods of 
examining candidates using an online database. 

The survey has so far been successful in finding 93 new 
pulsars, as well as 10 new RRATs and 1 RRAT candidate 
which were discovered using the algorithm described in 
this work. Among the newly discovered pulsars are 4 
millisecond pulsars which have been included in pulsar 
timing arr ays for the purpos e of detecting gravitational 
waves (e.g. Hobbs et al.|2010), a double neutron star sys¬ 
tem (Lynch et al., in prepj, and a binary system wit h an 


unusual degener ate companion (PSR J1816+4510; Ka- 
plan et aL]|2012 ). 


2.3. Survey sensitivity 

The sensitivities of the Drift-scan and the GBNCC sur¬ 
veys t o reg ular pulsars ar e prese nted by [Lynch et al.| 
(|2013|) and [Stovall et al.j ([2014[), respectively. Here 
we compute thes e surveys’ sensitivities to sin gle-pulse 
sources, following [Cordes fc McLaughhn[ ( [2003[ ) : 




/3(S/N)5(Tsys + Tsky) / IF, 


GWi 


ripAf ’ 


( 1 ) 


where Si is the intrinsic flux density of the pulse, P is 
a factor accounting for sensitivity losses due to digitiza¬ 
tion, (S/N)^ is the measured signal-to-noise ratio of the 
broadened pulse, Tgys is the system temperature at the 
frequency of observation, Tgky is the sky temperature, G 
is the telescope gain, Wi and IU 5 are the intrinsic and 
broadened pulse widths, respectively, Up is the number 
of summed polarizations, and A/ is the bandwidth. 

The detected pulse width, H4, depends not only on 
the intrinsic pulse width, but also on ISM propagation 
effects such as dispersion and scattering, which smear 
the pulse. Although dedispersion and a large number of 
frequency channels are used in order to mitigate pulse 
smearing due to dispersion, it is impossible to correct for 
the dispersion smearing across an individual frequency 
channel. The broadened pulse width is given by 




chan 


+ ^2 
' ^sc, 


( 2 ) 


where tgamp is the sampling time, tchan is the disper¬ 
sive smearing within each frequency channel, and tgcatt 
is the scatter broadening time related to the effect of 
multi-path sc attering of si gnals due to irregularities in 
the ISM (e.g. Rickett[[199Q[). An expre ssion for Ghan can 
be found in |Lorimer ki Kramer (2004); scaling to Drift- 
scan/ GBN(jC""dEservm^Trequency and channel band- 
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Fig. 1.— Sky coverage of the GET Drift-scan and GBNCC surveys, plotted in Equatorial coordinates. The gray strips indicate the 
areas of sky observed in the Drift-scan, while in green are the areas observed in the GBNGG survey. Discoveries from each survey are also 
indicated: pulsars are denoted by black circles, and RRATs are denoted by red stars. The sky coverage and discoveries of the GBNGG are 
current up to 2015 March. 


width, 


^chan — 4.7 flSX 


^/chan 

0.0244 MHz 


/ 


350 MHz 


-3 


DM 


An approximation for tgcatt is given by Cordes (2002), 


pc cm 

( 3 ) 


-3 


log(^^) =-3.59+ 0.129 log DM 

+ 1.02(logDM)2-4.41og(^^^ , (4) 

with A/chan the channel bandwidth and / the central 
observing frequency. 

For a minimum detected signal-to-noise ratio (S/N)^ 
of 6 and an assumed intrinsic pulse width, we can com¬ 
pute each survey’s sensitivity to single-pulse sources by 
substituting the survey-specific parameters from Table 
into Equation!^ We plot the resulting sensitivity as a 
function of DM in Figure for intrinsic pulse widths 
of 2 ms, 10 ms (the typical width of the RRAT pulses 
we detected), and 50 ms, for the Drift-scan and GBNCC 
surveys. 


3. RRAT SEARGH ALGORITHM 

The data from the Drift-scan and GBNCC surveys 
were analyzed using standard tools from the PRESTO 


suite of pulsar search software (Ransom 200l[p^ This 


analysis includes RFI removal, dedispersion, smgl e-pnlse 


searches, an d periodicity searches, a s described in[Lynch| 
et ah I ( 2013) for the Drift-scan and Stovall et al. poTIf 
lor CBNCC. RRATtrap, the single-pulse silting algorithrn 
developed here, was then applied to the outputs of the 
PRESTO single-pulse search routine in order to search for 
RRATs. We describe our algorithm and present its re¬ 
sults for the Drift-scan and GBNCC surveys. 


http://WWW.cV.nrao.edu/~ sransom/pre st o 



Fig. 2.— Minimum detectable flux density, ^min, of single pulses 
as a function of dispersion measure, for intrinsic pulse widths Wi 
of 2, 10, and 50 ms, plotted for the Drift-scan (blue) and GBNGG 
(green) surveys. The green band shows the scatter due to the 
variance in measured pulse broadening (Equation® for the Wi = 
50 ms GBNGG curve, with similar bands applicable for the other 
curves. 


3.1. Single-pulse sifting 

The data are searched for bright single pulses using 
presto’s single_pulse_search .py. This tool convolves 
the dedispersed and downsampled time series with box¬ 
car filters of various widths (namely 1, 2, 3, 4, 6, 9, 14, 
20, 30, 45, 70, 100, and 150 time bins) and records events 
which exceed the mean by an imposed threshold, her e 5 
standard deviations (see Cordes & McLaughlin|2003 for 
details on matched filtering). s ingle_pulse_search. py 
does not suffer from t he ‘root 2 problem’ reported by 
Keane & Petroff (2015), in which the sensitivity to pulses 
varies depending on the alignment of events with time 
series samples. This is an issue for programs which 
use power-of-two downsampling during the single-pulse 
search, rather than downsampling the data before search¬ 
ing for pulses and using matched filtering. 
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TABLE 1 

Observing modes used in RRAT discovery and follow-up observations. 


Parameters 

GBT Drift-scan 

GBNCC 

GBT 350 
follow-up 

GBT 820 
follow-up 

LOFAR 

Center freq. (MHz) 

350 

350 

350 

820 

150 

Bandwidth (MHz) 

50 

100 

100 

200 

80 

Number of freq. channels 

2048 

4096 

2048 

2048 

6400 

Sampling time {gs) 

81.92 

81.92 

81.92 

81.92 

655.36 

Obs. duration (min) 

2.3 

2 

10-15 

6-15 

10-30 

Beam FWHM (arcmin) 

36 

36 

36 

15 

5 

Backend 

Spigot 

GUPPI 

GUPPI 

GUPPI 

Blue Gene/P 

Analog-to-digital conversion factor, /3 

1.16 

1 

1 

1 

1 

System temperature (K)^ 

23 

23 

23 

23 

400^ 

Gain (K/Jy) 

2 

2 

2 

2 

5 

Sensitivity limit (mJy) 

260 

160 

160 

45 

80^ 


Note. — The sensitivity limit quoted above is for single-pulse sources, and is computed as in Section [2.3| 
for a DM of 100 pccm“^. 


The system temperature we report does not incl ude the sky tempe rature. 
The average receiver temperature is reporte d by 


Arts et al. 


The sensitivity for LOFAR is obtained from [Stappers et al. 


(20111 after scaling to the full available band¬ 
width, and assuming a 10-ms pulse width as in tRe calculations tor the GBT observations. 


(2013). 


A list and 4-5 plots of all S/N > 5 candidate single¬ 
pulse events in the beam are then saved, with each plot 
covering a portion of the DM range searched. Since 
single-pulse plots are produced for every pointing, the 
number of resulting plots for a full survey is very large 
(e.g. 4 X 30000 for the Drift-scan, and 5 x 66000 for GB- 
NCC so far), making their examination a tedious task 
with potentially subjective outcomes. We have devel¬ 
oped an automated sifting algorithm in order to identify 
RRAT candidates and astrophysical pulses in the output 
of single-pulse searches, eliminating the need for man¬ 
ual inspection of each diagnostic plot produced. This 
algorithm is described below and illustrated in Figure 
An example output plot from the algorithm is shown in 
Figure 


3 . 1 . 1 . Algorithm design 

The single-pulse sifting algorithm is designed to take 
advantage of properties that characterize astrophysical 
signals, and use these to distinguish astrophysical sig¬ 
nals from other signals. Our algorithm is based on the 
following concepts: 


(a) A bright signal will be detected over a range of DMs, 
with the strongest detection at the optimal DM and 
weaker detections above and below this DM due to 
smearing of the pulse when it is dedispersed at incor¬ 
rect DMs. The expected S/N fall-off with in correct 
DM is given by |Cordes & McLaughlinj (|2003|) , Equa¬ 
tion 12 . 


(b) Since signals are strongest at the optimal DM, we 
expect that signals of terrestrial origin (namely RFI) 
will peak at a DM of 0 pc cm“^. We can thus classify 
any signals that peak at DM ^ 0 pccm“^ as not 
astrophysical and reject them. 

(c) RFI that is only present in a narrow range of frequen¬ 
cies will not be subject to dispersion effects, and thus 
will appear consistently bright over a very large range 
of DMs. 

Concept 0 means that a given pulse, whether astro- 
physical or not, will be associated with many statisti¬ 


cally significant “single-pulse events” that will be found 
in the single-pulse search. These events will be spread 
over a small range of DMs, and will occur at approxi¬ 
mately the same time. The first step in our algorithm 
is thus to group events that belong to the same pulse 
by checking whether they satisfy this criterion, that is, 
lie within some small window of DM and time. Once 
all single-pulse events in a beam have been divided into 
groups, we examine each group’s collective properties in 
order to decide whether it behaves like an astrophysical 
pulse, and rate it based on these results. 

The first test we employ is related to group size. If a 
group has too few (< 20 , a tunable parameter) events, 
we classify it as noise, and assign it a rank of 1. We found 
that this minimum group size was sufficient in order to 
remove much of the noise, while still maintaining sensi¬ 
tivity to the weakest pulses = 6 ). Moreover, 

examination of detections of known pulses has shown us 
that this group size is a fairly conservative (i.e. easy for 
a real pulse to reach) criterion. However, a potential fu¬ 
ture improvement might be to relate the minimum group 
size for each pulse to its width and peak S/N, since these 
affect the number of events that will be associated with 
a given pulse. 

Next, we examine the S/N vs. DM behavior of the re¬ 
maining groups. From concept Q, we expect that a 
group of events that is caused by RFI will have a peak 
S/N at DM ^ 0 pccm“^. Thus, any group that satisfies 
this criterion, namely DM(S/N^^^) < DM^in, is classi¬ 
fied as RFI, which we define by a rank of 2. Here we use 
DMmin = 2 pccm“^. 

We then again make use of concept Q by looking for 
groups whose S/N peaks at some non-zero DM and de¬ 
creases above and below that DM. We divide each group 
into five intervals in DM (numbered 1-5 from low to 
high DM), with each interval containing one fifth of the 
group’s events. We refer to the maximum S/N of the 
events in interval i as S/N^. If the maximum S/N of the 
group occurs in one of the inner intervals {i = 2 ,3,4), and 
S/N^ drops off monotonically on either side of the peak 
(e.g., S/N^ < S/N 2 < S/N 3 and S/N 3 > S/N^ > S/N^), 
the group is classified as a “good” astrophysical pulse 
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and given a rank of 4. Moreover, if such a group’s peak 
S/N exceeds an imposed threshold, = 8 , the 

group is classified as “excellent” and given a rank of 5. 

If a group only somewhat follows this behavior but 
deviates from the expected trend (e.g. S/N 2 < S/N 3 and 
S/N 3 > S/N 4 , but S/N^ > S/N 2 or S/N 5 > S/NJ, it is 
given a rank of 3 and called “fair.” If instead a group has 
mostly constant S/N and spans a large range of DMs (in 
this case a range > 100 pccm“^, a tunable parameter), 
we conclude by concept Q that it is RFL Finally, any 
group with S/N behavior that does not fall into one of 
these categories is labelled “other” and given a rank of 0 . 

The values we used for the tunable parameters (max¬ 
imum DM span and minimum number of events per 
group) were chosen after testing various values on a sub¬ 
set of data, and were found to provide a good balance 
between the number of pulsars identified and the total 
number of candidates. These parameter values will de¬ 
pend on the observing frequency and RFI environment of 
a given pulsar survey. We find that on average, approxi¬ 
mately 63% of all pulses are categorized as RFI, 11% as 
fair, 13% as good, 0.1% as excellent, and 13% as other. 

For each beam, colorized single-pulse plots are pro¬ 
duced, with colors corresponding to group ratings (see 
Figure]^. A summary text file is also saved, containing 
each group’s characteristics including its time, minimum 
and maximum DM, pulse width, peak S/N, and assigned 
rating. These can in principle be used to sort candidate 
pulses, for example by assigning the characteristics dif¬ 
ferent weights or applying further decision algorithms. 

Once the beams have all undergone this sifting algo¬ 
rithm, those that have been flagged as having “excellent” 
pulses are visually examined. For those events that in¬ 
deed look astrophysical, we then generate frequency vs. 
time plots that show arrival times of the signal in dif¬ 
ferent frequency bins throughout the band. Since astro- 
physical signals experience a frequency-dependent ( 1 //^) 
dispersive delay while propagating through the interstel¬ 
lar medium, we expect to see this signature if the de¬ 
tected pulse is indeed astrophysical. We thus use these 
plots and our knowledge of the effects of the ISM as a 
final means of testing the astrophysical nature of a sig¬ 
nal, before deciding whether it is a candidate worthy of 
observational follow-up. 

3 . 1 . 2 . Algorithm performance 

We can assess the performance of our algorithm by 
quantifying its false positive and false negative rates. 
Moreover, we can examine the algorithm’s usefulness by 
comparing the number of beams that it flags as poten¬ 
tially interesting to the total number of beams analyzed, 
since the latter is the number of beams that would have 
to be visually inspected if the sifting algorithm were not 
used. We find that on average, the algorithm identifies 
about 10 % of all beams as containing “excellent” astro- 
physical pulses. This means that we only have to exam¬ 
ine 10 % of all diagnostic single-pulse plots, which is a 
significant improvement. 

The false positive rate may be computed based on 
the number of beams identified as containing “excellent” 
pulses which, when inspected visually, are not found to 
contain astrophysical pulses. This misidentification is 
likely due to noise and unmodelled behavior of RFI that 
leads to clusters of single-pulse events that satisfy the 


imposed criteria for astrophysical pulses. We find that 
this rate is approximately 90%, meaning that of every 10 
beams identified, one will contain a pulsar (either known 
or new). If the algorithm were not used and the diagnos¬ 
tic single-pulse plots for every beam were examined, we 
would have to look through 100 beams in order to find 
this one pulsar. Nonetheless, we are currently working 
on improving the false positive rate in order to make the 
algorithm more effective, and further reduce the number 
of diagnostic plots that must be examined. 

Since we use our algorithm to identify which beams 
should be inspected, it is also possible that some pul¬ 
sars are missed. Quantifying this false negative rate is 
more difficult. Ideally, we would compute this rate by 
examining the output of the code for all beams contain¬ 
ing known RRATs and single-pulse visible pulsars, or by 
performing a blind injection of synthetic RRAT signals 
into a sample of data files. Instead, here, we obtain a 
rough approximation of this rate by examining a sample 
of beams which were found to contain RRATs and pul¬ 
sars (both new and known), and check the output of the 
code to see how many of their pulses were missed and 
not rated as “excellent.” We find that this rate is ^20%. 
If we loosen our criteria and count pulses which were 
not rated as “excellent” or “good,” this rate decreases 
to ^ 10 %, though we note that in practice only beams 
with top rated (“excellent”) groups were examined in 
this work. This rate implies that there is a 20% chance 
of missing an individual pulse. Note, however, that this 
value is an average, and in practice depends strongly on 
properties such as pulse shape and width, DM, and S/N. 

The false negatives in our analysis were most com¬ 
monly caused by S/N vs. DM curves that did not con¬ 
form to expectations. Qualitatively, we expect the S/N 
vs. DM curve to have a single maximum at the optimal 
DM and decrease monotonically with increasing DM dif¬ 
ferences (Cordes & McLaughlin||2003). However, this is 
based on the assumption that the pulse shape has a sin¬ 
gle component. In practice, complex multi-component 
pulses produce asymmetric or multi-peaked S/N vs. DM 
curves (see, e.g., the ‘Fair’ panel in Figurethat, if suf¬ 
ficiently extreme, result in false negatives. Another case 
in which our S/N vs. DM curve test may fail is for weak 
pulses near the noise threshold. In this case, the fluctu¬ 
ations in S/N due to noise are significant compared to 
the peak S/N of the pulse. We note, however, that both 
categories of false negatives described here still contained 
many “fair” pulses. While we can make our criteria less 
stringent (e.g. by looking at all “fair”, “good” and “ex¬ 
cellent” pulses) in order to decrease the false negative 
rate, doing so currently would immensely increase the 
false positive rate, thus lowering the effectiveness of our 
algorithm. However, we are working toward improving 
the algorithm in order to decrease the false positive and 
negative rates, which will also allow us to probe weaker 
pulse candidates. 


3 . 2 . RRAT discoveries 

We applied our RRAT search algorithm to the Drift- 
scan survey, and to 21300 beams of the GBNCC survey, 
representing about a third of the GBNCC data taken 
thus far. This resulted in the discovery of 21 new RRATs 
and 6 RRAT candidates. The positions of these RRATs 
are indicated on the survey sky coverage map, in Fig- 
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Fig. 3.— RRATtrap algorithm flowchart, illustrating the rating process applied to candidate single-pulse events, as described in Sect ion [3.1| 
The threshold S/N is plotted as a dashed line in the three right plots for illustration. 


ure^ We present the properties of each RRAT in Ta¬ 
ble 0 including position, period, DM, and the survey in 
whi^ it was found. We note that these parameters were 
obtained through follow-up observations as will be de¬ 
scribed in Section and the refined values are the ones 
reported. The RRAT candidates, as well as the number 
and duration of re-observations for each one, are listed 
in Table The discovery plots for the candidates are 
shown in Figure We note that these candidates are 
all fairly close in DM; some of them may therefore be 
caused by RFI, or may simply be coincidentally close in 
DM. We define a candidate as any source that was only 
detected in the initial survey observation, and was not 
re-detected in any observations thereafter. These could 
thus be spurious astrophysical signals, or they could be 
genuine RRATs with low burst rates. 


4. RRAT FOLLOW-UP OBSERVATIONS AND ANALYSIS 

We conducted follow-up observations of our RRAT 
discoveries between February 2012 and January 2015 
using the GET and LOFAR telescopes. Observations 
with the GET were conducted at center frequencies of 
350 MHz and 820 MHz, with bandwidths of 100 MHz 
and 200 MHz, respectively. Integration times ranged 
from 6-15 minutes. Observations with LOFAR used a 
center frequency of 150 MHz, a bandwidth of 80 MHz, 
and integration times of 10-30 minutes. These and other 
observation parameters are summarized in Table for all 
observing modes used. These observations allowed us to 
refine positions, periods, DMs, and burst rates for our 
RRATs, using the methods described below. We also ob¬ 
tained phase-coherent timing solutions for four RRATs. 















































TABLE 2 

New RRAT Discoveries in Green Bank Telescope 350 MHz Surveys. 


RRAT 

R.A. 

(hh:mm:ss) 

Dec. 

(dd:mm) 

P 

(s) 

DM 

(pc cm“^) 

Survey 

J0054+69 

00:54:28(25) 

+69:26(2) 


90.3(2) 

GBNGG 

J0103+54 

01:03:37(20) 

+54:02(5) 

0.354304(6) 

55.605(4) 

GBNGG 

J0201+7005* 

02:01:41.344(7) 

+70:05:18.11(6) 

1.349184471847(9) 

21.029(2) 

GBNGG 

J0332+79 

03:32:45(20) 

+79:10(5) 

2.05621(6) 

16.67(2) 

GBNGG 

J0447-04 

04:47(1) 

-04:35(18) 

2.18819(2) 

29.83(4) 

Drift 

J0545-03 

05:45(2) 

-03:10(36) 

1.07393(2) 

67.2(4) 

Drift 

J0957-06 

09:57(2) 

-06:17(36) 

1.72370(8) 

26.95(2) 

Drift 

J1126-27 

11:26(2) 

-27:37(36) 

0.358161(8) 

26.860(7) 

GBNGG 

J1153-21 

11:53(2) 

-21:18(36) 

2.34348(5) 

34.8(1) 

GBNGG 

J1332-03 

13:32(2) 

-03:26(36) 

1.10640(5) 

27.1(2) 

Drift 

J1439+76 

14:39(2) 

+76:55(36) 

0.947903(2) 

22.29(2) 

GBNGG 

J1538+2345* 

15:38:06.07(2) 

+23:45:04.0(2) 

3.44938495332(9) 

14.909(1) 

Drift 

J1611-01 

16:11(2) 

-01:28(36) 

1.29687(2) 

27.21(7) 

Drift 

J1705-04 

17:05(2) 

-04:41(36) 

0.23748(2) 

42.951(9) 

Drift 

J1915-11 

19:15(2) 

-11:30(36) 

2.1770(2) 

91.06(8) 

Drift 

J1944-10 

19:44(2) 

-10:17(36) 

0.409135(1) 

31.01(3) 

Drift 

J1956-28 

19:56(2) 

-27:53(36) 

0.2600144(6) 

45.69(1) 

GBNGG 

J2007+20 

20:07(2) 

+20:21(36) 

4.634(7) 

67.0(4) 

Drift 

J2105+6223* 

21:05:12.93(2) 

+62:23:05.5(1) 

2.30487883766(4) 

50.75(8) 

GBNGG 

J2311+67 

23:11(1) 

+67:05(18) 

1.9447(2) 

97.1(2) 

GBNGG 

J2325-0530* 

23:25:15.3(1) 

-05:30:39(4) 

0.868735115026(9) 

14.966(7) 

Drift 


Note. — Periods that could not be reliably constrained are not reported (see text). Uncer¬ 
tainties on measurements are given in brackets, representing the uncertainty on the last reported 
digit. Some source positions are bet ter d etermined than others thanks to position refinement 
observations j as described in Section |4.1| and others are even better determined using timing 
(see Section |4.4| ) . Sources for which we nave full timing solutions are indicated by an asterisk 
and are presented in Table ^ 


TABLE 3 

Unconfirmed RRAT Candidates in Green Bank Telescope 350 MHz Surveys. 


RRAT 

R.A. 

(hh:mm) 

Dec. 

(dd:mm) 

p 

(s) 

DM 

(pccm“^) 

Number of pulses 

Peak fiux density 

(Jy) 

Survey 

Re-observations 

(min) 

J0441-04 

04:41 

-04:18 


20.0 

1 

0.4 

Drift 

15 + 3 + 15 + 10 

J0513-04 

05:13 

-04:18 


18.5 

1 

1.0 

Drift 

10 + 10 + 10 

J0614-03 

06:15 

-03:29 

0.136 

17.9 

3 

0.9 

Drift 

6 + 10 + 10 

J1059-01 

10:59 

-01:02 


18.7 

1 

0.8 

Drift 

10 + 10.5 

J1336-20 

13:36 

-20:34 

0.184 

19.3 

3 

0.8 

Drift 

10 + 10 + 10 

J1354+25 

13:54 

+24:54 


20.0 

1 

0.1 

GBNGG 

10 


Note. — True period values may be integer fractions of the tabulated values, as described in Section |4.3| Periods for 
sources that showed two or fewer pulses per observation could not be reliably constrained and are not reported. No detailed 
period and DM refinement, nor a robust determination of uncertainties, was done for these sources since they are pending 
confirmation. The peak fiux density of each candidate was computed using the observed pulse width (assumed to be the 
optimal box-car width used in the search), maximum S/N, and telescope parameters. A summation of integration times in 
the ‘Re-observations’ column indicates separate observation epochs. 


We describe this timing analysis below and present the 
results in Table [U 

LOFAR is made up of many dipole antennas which 
are arranged hierarchically in stations, with core sta¬ 
tions densely packed near Exloo, the Netherlands, re¬ 
mote stations spread throughout the Netherlands, and 
international stations in Germany, France, Sweden, and 
the United Kingdom. For our LOFAR observations, we 
used all core high-band antenna (HBA; 110-190 MHz) 
stations. The LOFAR core stations subtend a maximum 
baseline of approximately 2 km and create a tied-array 
beam with a FWHM of approximately 5' (at 150 MHz). 
Typically 20-23 core stations were used in each LOFAR 
observation. When observing using multiple stations, 
one can take advantage of multi-beaming in order to 
form and observe several beams on the sky simultane¬ 


ously. Moreover, adding signals from stations coherently 
results in a tied-array beam, which has high spatial res¬ 
olution and increased sensitivity. We used these configu¬ 
rations for our observations with LOFAR, as we describe 
below. More details on available observing modes with 
LOFAR are presented in van Haarlem et al. (2013) and 


Stappers et al. (2011). 


4.1. Position refinement 

Positions were refined using standard gridding with the 
GET at 820 MHz for some sources, and using LOFAR 
beamforming for other sources, reducing positional un¬ 
certainties to 15' (GET) or 5' (LOFAR). In particular, 
LOFAR was especially useful for the gridding observa¬ 
tions, since we were able to use all core high-band an¬ 
tenna stations to form 61 simultaneous tied-array beams. 
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Single Pulse Results for J1537+2350 
RA: 15:37:47.3637 Dec: 23:50:51.1319 
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Fig. 4.— Example of RRATtrap single-pulse sifting output, for 
newly discovered RRAT J1538+2345. Each candidate single-pulse 
event is plotted, with the point size proportional to the event’s 
S/N. Radio frequency interference is identified in red, whereas 3 
pulses are successfully identified in purple, and another weaker 
pulse in blue. These pulses all appear strongest at a common DM of 
~15 pccm“^. Some RFI is present but not successfully identified 
due to its irregular shape (seen in black and blue, at t = 20 — 30 s). 
Furthermore, a fifth pulse is visible by eye, at the same DM as 
the identified pulses and at t = 38 s, but is not identified by the 
code because its S/N does not follow the prescribed behavior (see 
Figure]^. 

each of width 5', and arranged in four concentric rings, 
thus efficiently covering the 350-MHz GBT beam. This 
setup is particularly useful for localizing RRATs, com¬ 
pared to sequential gridding observations which require a 
longer integration time overall and are made challenging 
by RRATs’ sporadic emission and pulse amplitude vari¬ 
ations. Figureshows an example of the LOFAR config¬ 
uration for gridding observations of one of the RRATs. 

In the analysis of LOFAR gridding observations, each 
beam is processed independently using standard pul¬ 
sar analysis tools, and summary plots are then pro¬ 
duced for the observation to present results for all tied- 
array beams. These diagnostic plots include detec¬ 
tion ‘heatmaps’ of two types: The periodicity heat map 
shows the relative detection significance when each beam 
is folded at the nominal period of the RRAT, and 
the RRAT heatmap compares the single-pulse detection 
strength at the DM of the RRAT in each beam. An ex¬ 
ample of the single-pulse heatmap for the observation of 
RRAT J0054+69 is seen in Figure and the correspond¬ 
ing detection is shown in Figure |R The RRAT heatmap 
provides an overview of the single-pulse results without 
requiring the inspection of a single-pulse plot for every 
beam, though the latter may still be needed when the 
detection is weak or when RFI is present. 

4.2. Burst rate calculations 

An important characteristic of RRATs, and one of the 
properties that determines whether or not a source can 
be timed, is the source burst rate. For every observation 
of a given source, we compute the burst rate by divid¬ 
ing the number of detected pulses by the observation 
duration. Note that observed burst rates are observing- 
setup-dependent, as a telescope’s sensitivity will affect its 


ability to detect weaker pulses, and single pulses from a 
given source vary in amplitude (e.g. Burke-Spolaor et al. 
2012). RRATs may also behave difterently at difterent 
frequencies, depending on their spectrum. It is gener¬ 
ally expected that RRATs, like pulsars, are broadband 
sources, and emit over the entire radio spectrum. How¬ 
ever, if one finds significantly different burst rates at dif¬ 
ferent frequencies, some inferences regarding the source 
spectrum can be made. 

Moreover, since RRATs are erratic emitters, their mea¬ 
sured burst rates may change drastically between obser¬ 
vations using the same setup. It is thus important to use 
data from multiple, sufficiently long observations, and 
report the average burst rate, computed by dividing the 
number of detected pulses over all observations by the 
summed observation time. We estimate the uncertainty 
on the number of detected pulses using Poisson statistics. 

A summary of source detections and burst rates with 
each observing setup is presented in Table Since 
LOFAR’s sensitivity is elevation-dependent, the LOFAR 
burst rates of different RRATs are not directly compa¬ 
rable to one another. Moreover, many RRATs were not 
detected with LOFAR, which may be due to reduced 
sensitivity at low elevations (those sources with low de¬ 
clinations only reach a maximum elevation of ^30° at 
LOFAR), higher sky temperatures close to the Galactic 
Plane, or their spectra. 

4.3. Timing and period refinement 

When RRATs are initially discovered, their periods can 
be constrained by finding the greatest common divisor of 
the time differences between detected pulses. This can 
typically be done as long as three or more pulses are ob¬ 
served. The obtained value may then be an integer mul¬ 
tiple of the true period, but as more pulses are observed 
we can establish a period value with higher confidence. 
These period values are then refined through follow-up 
observations and timing techniques, as described below. 
Furthermore, if a sufficiently long data set is available, 
we can derive a phase-coherent timing solution. 

We have timed four RRATs, and refined the single- 
pulse-inferred period values for the other sources. For 
one RRAT, J0054+69, we observed an abundance of 
pulses, however we were unable to fit a periodicity to 
these pulses. We suspect that this may be due to a 
complex pulse profile, with detected single pulses cor¬ 
responding to different profile components. In fact, even 
for single-peaked average profiles, individual pulses from 
a RRAT are narrower and distributed across the enve¬ 
lope of the average pulse profile, introducing jitter which 
makes timing difficult. 

We start by extracting all single pulses from the data. 
The output of single_pulse_search.py provides ap¬ 
proximate times of detected pulses. One way of proceed¬ 
ing, then, for a source whose period is roughly known, 
is to divide the observation into time segments corre¬ 
sponding to single rotations (or ‘pulse archives’), and ex¬ 
amine those data segments near the approximate pulse 
times by eye until identifying bright puls es. We perform 
the data cho pping using the DSPSR tool (van Straten & 
Bailes||2011 ), which also allows to simultaneously dedis- 
perse the data ac cording to the RRA T’s DM. We then use 
PSRCHIVE tools (jHotan et al.||2004|) to make frequency- 
versus-phase plots for each pulse archive, where pulses 
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Single pulse results for 'GBT350drifL54238_0441-0418' 


Single pulse results for 'GBT350drifL54238_0513-0417' 


Source: unknown 
Telescope; GBT 
Instrument: SPIGOT 
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DEC (J2000): -04:18:13.3422 Sampling time: 163.84 ;t4S 
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Instrument: SPIGOT 
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Single pulse results for 'GBT350drifL54224_0004_0078' 

Source: unknown RA (J2000): 06:14:37.8153 N samples: 864000 

Telescope: GBT DEC (J2000): -03:28:44.9314 Sampling time: 163.84 us 
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Single pulse results for 'GBT350dnft_54274_1059-0102' 


Source: unknown 
Telescope; GBT 
Instrument: SPIGOT 
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Single pulse results for 'GBT350dnft_54294_1 336-2034' 


Single pulse results for *guppL56460_GBNCC43999_0022_000r 


Source: unknown 
Telescope; GBT 
Instrument: SPIGOT 
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Fig. 5. — Discovery plots of the six unconfirmed RRAT candidates listed in TableJ^ From left to right and top to bottom, these are 
candidates J0441—04, J0513—04, J0614—03, J1059—01, J1336—20, and J1354+25. Each plot is made up of four subplots, showing the 
number of events vs. S/N, the number of events vs. DM, the S/N vs. DM, and DM vs. time of detected events, with the plotted circle 
size corresponding to each event’s S/N. 
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Fig. 6. — Single-pulse detection heatmap for LOFAR gridding ob¬ 
servation of RRAT J0054+69, taken December 20, 2013, as part of 
our follow-up campaign. The heatmap compares detection signifi¬ 
cance of single-pulse events at the RRAT’s DM across all beams. 
The S/N values are summed for all events at the given DM, for 
each beam. A detection is evident in beam 3. The apparent de¬ 
tections in some of the other beams (e.g. 11, 40, 60) are due to 
RFI signals. Note that the circle sizes do not represent the size of 
the tied-array beams on the sky. The synthesized beam is shown 
in the bottom-left corner. 
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Fig. 7.— LOFAR detection of single pulses from RRAT 
J0054+69. Approximately 7 pulses at a DM of 90 pccm“^ are 
detected during a gridding observation, allowing the refinement of 
the position of this source to that of beam 3 of the 61 tied-array 
beam grid, thus lowering the position uncertainty from 36' to 5'. 
A bright RFI signal is seen at t ~ 900 s. 


can be identified. 

Once all single pulses for a given observation have been 
identified, we align the pulses by cross-correlation and 
add them together using psradd to create a high S/N 
pulse. We then use paas to interactively fit an analytic 
pulse profile template to the combined pulse. We use this 
as an initial template for the remaining observations (if 
available), but further refine the template once a timing 


solution is obtained (see Section 4.4). A time-of-arrival 
(TOA) is then computed for each single pulse, using 
the standard cross-correlation algorithm in the pat tool 
( [Taylor [l992 ). After extracting TOAs for each source, 
we use TEMPOand standard pulsar timing techniques 
to refine the period within a single observation for each 
source. 

We used these techniques to refine the periods of the 
14 RRATs for which we did not have many observations. 
We used a similar approach but different software for 
the four RRATs for which we had a t imin g baseline of 
observations, as we describe in Section [44 


4.4. RRATs with timing solutions 

We obtained phase-coherent timing solutions for four 
RRATs using similar techniques to those described 
above. For these RRATs, we had more data and pulse 
detections, and were thus able to obtain stable pulse 
profiles using high S/N observations or multiple obser¬ 
vations summed together. We therefore did not need the 
fine-grain control offered by DSPSR and PSRCHIVE, and 
instead used PRESTO tools for this analysis. We folded 
the data with prepfold while specifying a large num¬ 
ber of subintegrations, each corresponding to a few sec¬ 
onds, generated TOAs using get_T0As.py, and ignored 
all subintegrations that did not contain pulses. We found 
that two of the four RRATs were consistently detectable 
from folding the data, but this method worked well for all 
four nonetheless. We then used TEMP02 for the timing. 

The timing solutions for these RRATs are presented 
in Table [^ and their residuals are shown in Figure [^ 
We also plot these RRATs on the P — P diagram in 
Figure [^ Comparing these four RRATs to the previously 
known RRATs, we see that they lie roughly in the same 
region, with periods of order 1 s and magnetic fields of 
order 10^^ G. 

We did not have sufficient data in order to obtain tim¬ 
ing solutions for the other RRATs, and in particular were 
not able to detect many of them with LOFAR (as seen in 
Table 1^. Given sufficient telescope time, we should be 
able to obtain timing solutions for the remaining RRATs, 
for example using the GET. It is also possible that some 
of these sources would be detected in longer observations 
with LOFAR. There are also upcoming telescopes, such 
as the Ganadi an Hydrogen Intensit y Mapping Experi¬ 
ment (GHIME; [Bandura et al.[[2014[), which will provide 
daily coverage of the sky and will thus be very useful for 
observing and timing RRATs. 

4.5. DM determination 

The DM of a RRAT can be approximated by exam¬ 
ining the S/N behavior of a group of single-pulse events 
asso ciated with an individual pulse. As described in Sec¬ 
tion 3.1.1[ the S/N of a pulse will peak when the data are 
dedispersed at the optimal DM. We can thus estimate the 
DM as the value at which the RRAT’s pulses are bright¬ 
est. This method requires no additional computation as 
it uses data products produced by the data reduction 
pipeline. However, because the single-pulse search in the 
pipeline uses computationally efficient box-car filters as 
an approximation to the actual spiky shape of astrophys- 
ical pulses, the DM values determined from this method 


http://WWW.atnf.csiro.au/research/pulsar/tempo2 
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have some slight systematic errors that are difficult to 
quantify. 

In order to refine the DM of a RRAT, we first iden¬ 
tify all single pul ses using DSPSR and PSRCHIVE as de¬ 
scribed in Section [4~3l but without dedispersing the data. 
We then choose the brightest pulse, or a sum of several 
pulses, divide the frequency band into 4-6 subbands, and 
extract a TOA for each subband. Due to the disper¬ 
sion sweep, we expect that these TOAs will obey the 
1 //^ plasma dispersion law. We then input these TOAs 
into TEMP02 and fit for the DM, yielding a refined DM 
value. We fit this DM value using one observation for 
each source. 

5. DISCUSSION AND POPULATION COMPARISONS 

We now consider our new RRATs in the context of the 
previously known RRAT population, as well as the gen¬ 
eral radio pulsar population. We examine the observed 
burst rates for our RRATs, and discuss their implica¬ 
tions on survey detections. We then compare properties 
such as period and DM for RRATs and pulsars, in order 
to assert whether the underlying distributions are sta¬ 
tistically different. We consider the relevant detection 
selection biases, and describe how to conduct fair statis¬ 
tical comparisons despite them. Finally, we discuss the 
fact that we did not find any FRBs with relation to the 
published FRB rates. 

5.1. Burst rates 

As is evident from the new RRATs’ burst rates, which 
are recorded in Table and plotted in Figure \To\ there 
is a distribution of rates with values varying drastically 
from source to source. In fact, even for a given source, the 
measured burst rate may vary between observations due 
to small-number (Poisson) statistics, and the expected 
rate also varies with the observation’s sensitivity (includ¬ 
ing factors such as observing setup, source elevation, and 
RFI occupancy). It is therefore difficult to isolate the in¬ 
trinsic burst rate from the observed rate, and to know 
whether we have found the true average burst rate. The 
reported values should therefore be considered approxi¬ 
mate. 

These modest burst rates demonstrate the potential 
of missing a RRAT in a short observation, and suggest 
that the number of RRATs detected in a given surveyed 
region is only a fraction of the number therein. By ap¬ 
proximating RRAT pulse emission as a random process 
that follows Poisson statistics, we can estimate the true 
number of RRATs present within a certain region given 
the number of detected RRATs and an average burst rate 
for the population. 

The Poisson probability function is given by 

e~^ 

Pr(X = fc) = ^j-, (5) 

where, in this case, X is the number of detected pulses, 
Pr(X = k) is the probability of detecting k pulses, and 
A is the expected value of X given by the product of 
observation duration and mean burst rate. The proba¬ 
bility functions for different mean burst rates are plotted 
in Figure El where we have assumed a typical 2-min 
survey observation. For each burst rate, the probabil¬ 
ity of missing a RRAT over an observation is given by 


Ft{X = 0), whereas the probability of detecting one or 
more pulses is given by the sum Pr(X = /c), or 

simply I — Ft{X = 0). 

We can thus assume some mean burst rate for the 
RRAT population, and use the above to compute the 
probability of non-detection. In this way, we obtain an 
estimate of the fraction of the RRAT population missed 
in a given pulsar survey due to short observations and 
infrequent emission. It should be noted however that 
the observed burst rate distribution for our sources is in¬ 
herently biased, since these sources had sufficiently high 
burst rates to be detected. In fact, comparing the follow¬ 
up to discovery burst rate of each source, we see that 
the follow-up rates are generally lower, hinting that the 
survey burst rates are biased toward high values. We 
therefore cannot assume the average burst rate from our 
discoveries as that of the entire RRAT population. In¬ 
stead, we only illustrate this using several RRAT burst 
rates. 

With the observed burst rates tabulated in Table _ 
we first choose a rate in the low range of those va! 
ues, 30 pulses hr“^. We then find that for a 2-min 
observation, the probability of detecting zero pulses is 
about 37%. This implies that for every RRAT present 
in the surveyed area, there is only a 63% chance of de¬ 
tection, and we may thus apply a correction factor of 
1/0.63 = 1.6 to the number of detected RRATs to es¬ 
timate the true number of RRATs in the surveyed re¬ 
gion. This conclusion is of course highly dependent on 
the chosen mean burst rate. On the other hand, we can 
use the average observed burst rate from our sources, 
which is 86 pulses hr“^ (from follow-up observations). 
In this case, the probability of detecting zero pulses in 
a 2-min observation is 5.7%, yielding a correction factor 
of 1.06 to the expected number of RRATs in the region. 
This can thus be taken as a lower limit on the true cor¬ 
rection factor that should be applied. Another extreme 
for the assumed RRAT burst rate can be obtained from 
the average bur st rate of the Parke s Multi-beam Pulsar 
Survey (PMPS; |Keane et al.||20I0|) RRATs, which had 
longer observations and whose rates were thus better de¬ 
termined. From the RRATalog, this average burst rate 
is 4.2 pulses hr“^. Using this rate in the above analysis, 
we find that the probability of non-detection in a 2-min 
observation is 87%, yielding a correction factor of 7.7 to 
the number of RRATs in our survey region. 

We can also compute the minimum burst rate required 
for a RRAT to be detected with a given probability. 
For example, if we set the probability of detection to 
99%, we find that the minimum burst rate required is 
138 pulses hr“^ for a 2-min observation. This value is in¬ 
dicated on the burst rate distribution of our new RRATs 
in Figure We note that there are many detected 
sources below this cutoff, which is statistically unlikely 
and suggests that there is a large population of sources 
with low burst rates yet to be discovered in our survey 
region. 


5.2. Periods 

We now examine potential biases influencing the ob¬ 
served spin period distributions of detected pulsars and 
RRATs. The first point to address is the potential ambi¬ 
guity in measured periods. While the period of a pulsar 
found in a periodicity search is easily established, for 
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Fig. 8.— Timing residuals for the four RRATs. Residuals from GET 820-MHz data are plotted in blue, whereas those from LOFAR 
observations are plotted in red. Note the different scale for each subplot. 


TABLE 4 

Timing parameters for four RRATs. 



RRAT J0201+7005 

RRAT J1538+2345 

RRAT J2105+6223 

RRAT J2325-0530 

Measured parameters 

Right ascension (J2000) 

02:01:41.344(7) 

15:38:06.07(2) 

21:05:12.93(2) 

23:25:15.3(1) 

Declination (J2000) 

+70:05:18.11(6) 

+23:45:04.0(2) 

+62:23:05.5(1) 

-05:30:39(4) 

Spin frequency (Hz) 

0.741188488948(5) 

0.289906755417(8) 

0.433862285366(7) 

1.15109885937(1) 

Spin frequency derivative (Hz/s) 

-3.0294(8) X 10-1^ 

-5.79(1) X 10-1® 

-9.824(1) X 10-1® 

-1.363(2) X 10-1® 

Dispersion measure (pccm“^) 

21.029(2) 

14.909(1) 

50.75(8) 

14.966(7) 

Derived parameters 

Galactic longitude, 1 (°) 

128.89 

37.32 

99.79 

75.58 

Galactic latitude, b (°) 

8.03 

52.39 

10.18 

-60.20 

Spin period, P (s) 

1.349184471847(9) 

3.44938495332(9) 

2.30487883766(4) 

0.868735115026(9) 

Spin-down rate, P (s/s) 

5.514(1) X 10-1^ 

6.89(1) X 10-1® 

5.219(6) X 10-1® 

1.029(2) X 10-1® 

Magnetic field, B (G) 

2.8 X 1012 

4.9 X 1012 

3.5 X 1012 

9.6 X loll 

Spin-down luminosity, E (erg s“^) 

8.9 X 10^1 

6.6 X 10®° 

1.7 X 10®i 

6.2 X 10®i 

Gharacteristic age, Tc (yr) 

3.9 X 10® 

7.9 X 10® 

7.0 X 10® 

1.3 X 10^ 

DM distance, D (kpc)^ 

1.1 

1.2 

2.6 

0.7 

Observation parameters 

Discovery (MJD) 

55163 

54243 

55201 

54240 

Timing epoch (MJD) 

56777 

56745 

56774 

56774 

Start epoch (MJD) 

56521 

56455 

56514 

56514 

Finish epoch (MJD) 

57035 

57036 

57036 

57036 

Number of TO As 

204 

213 

16 

132 

RMS post-fit residuals {fis) 

610 

3203 

523 

1165 


Note. — ^ DM distances are as implied by the |Cordes &: Lazio] ( |200^ NE2001 model and have typical uncertainties of 25%. 
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Fig. 9.— P — P diagram for all known neutron stars outside of globular clusters. The four new RRATs are indicated by stars, the 
previously known RRATs with measured P by large circles, magnetars by crosses, XDINS by X’s, and other neutron stars by dots. Lines 
of constant magnetic field are dashed and lines of constant characteristic age are dotted. The solid line represents the ‘death line’ (using 
Equation 9 of|Chen & Ruderman|1993|), beyond which we do not expect to find pulsars, but note that the exact location of the death line is 
uncertain since tJie emis sion mechanism of pulsars is not fully understood. Values for pulsars were t aken from the ATNF Pulsar Catalogue 
( [Manchester et al.|2Q05| >, for magnetars from the McGill Magnetar Catalog ( jOlausen fc Kaspi|2014| >, and for RRATs from the RRATalo^ 


^ http://astro.phys.wvu.edu/rratalog 


TABLE 5 

RRAT BURST RATES FROM DISCOVERY AND FOLLOW-UP OBSERVATIONS WITH THE GBT AND LOFAR TELESCOPES. 


RRAT 

Discovery 
Burst rate (hr“^) 

t (min) 

GBT 350 
Burst rate (hr“^) 

t (min) 

GBT 820 
Burst rate (hr“^) 

t (min) 

LOFAR 
Burst rate (hr“^) 

t (min) 

J0054+69 

300 ± 95 

2 

186 ± 33 

10 

285 ± 21 

40 

6 ± 2 

70 

J0103+54 

390 ± 108 

2 

60 ± 15 

15 

0 

30 

3 ± 2 

45 

J0201+7005 

180 ± 73 

2 

24 ± 10 

15 

256 ± 22 

30 

63 ± 7 

75 

J0332+79 

30 ± 30 

2 

65 ± 12 

25 

66 ± 10 

40 

14 ± 3 

70 

J0447-04 

103 ± 51 

2 

41 ± 10 

25 

78 ± 22 

10 

0 

15 

J0545-03 

77 ± 45 

2 

42 ± 16 

10 

- 

- 

0 

15 

J0957-06 

180 ± 68 

2 

138 ± 29 

10 

- 

- 

0 

10 

J1126-27 

180 ± 73 

2 

48 ± 17 

10 

- 

- 

- 

- 

J1153-21 

150 ±67 

2 

60 ± 24 

10 

- 

- 

- 

- 

J1332-03 

51 ± 36 

2 

20 ± 5 

40 

- 

- 

- 

- 

J1439+76 

450 ± 116 

2 

406 ± 48 

11 

0 

7 

0 

15 

J1538+2345 

129 ± 58 

2 

77 ± 14 

23 

- 

- 

66 ± 7 

75 

J1611-01 

51 ± 36 

2 

60 ± 19 

10 

- 

- 

0 

15 

J1705-04 

26 ± 26 

2 

32 ± 11 

15 

- 

- 

0 

15 

J1915-11 

26 ± 26 

2 

108 ± 25 

10 

- 

- 

0 

15 

J1944-10 

180 ± 68 

2 

78 ± 13 

27 

- 

- 

0 

15 

J1956-28 

120 ± 60 

2 

66 ± 20 

10 

- 

- 

- 

- 

J2007+20 

77 ± 44 

2 

11 ± 4 

34 

- 

- 

0 

15 

J2105+6223 

30 ± 30 

2 

- 

- 

- 

- 

12 ± 3 

75 

J2311+67 

60 ± 42 

2 

- 

- 

168 ± 32 

10 

0 

75 

J2325-0530 

103 ± 51 

2 

46 ± 9 

33 



52 ± 8 

45 


Note. — A field containing indicates the source was not observed for follow-up with that setup, whereas 0 indicates the source 
was observed and not detected. Average source burst rates and uncertainties are computed as described in the text. Not e that the 
observed burst rate for a given source may vary significantly between different observing setups, as described in Section |T^ 
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Fig. 10. — Burst rate distribution of our new RRATs. The rates 
plotted are those inferred from the GBT 350-MHz discovery and 
follow-up observations, corresponding to columns 2 and 3 of Ta- 
ble^ respectively. The 99% Poisson burst rate cutoff is also shown, 
indicating the rate below which it is unlikely to detect RRATs in 
a 2-min observation. 
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Fig. 11.— Poisson probability function showing the probability 
of pulse detection over a 2-min observation, for different assumed 
RRAT burst rates. 


RRATs there is initially some ambiguity in the inferred 
period due to the way in which it is computed. That 
is, one may argue that the computed period is an inte¬ 
ger multiple of the true spin period, and that the long 
periods observed for RRATs are due to this incorrect 
measurement. However, as a RRAT is monitored over 
time, the probability of this occurring diminishes, since 
a period can be established with high confidence as more 
pulses are detected. 

One might also expect that a RRAT with a short pe¬ 
riod will have more chances of detection over a given 
observation compared to a slower source, since it com¬ 
pletes more rotations in the same amount of time. How¬ 
ever, it is also possible that this source will instead be 
detected as a regular pulsar through a periodicity search 
if its burst rate is sufficiently high. For example, some 
of the RRATs found in our 2-minute survey integrations 
may be detected in the Fast Fourier Transforms (FFTs) 
of longer observations, as was sometimes the case in our 


^ 20-minute follow-up observations. The classification of 
a source as a RRAT or normal pulsar is therefore related 
to the number of pulses detected, which i s a function 
of period, burst rate, and integration time. |McLanghlin| 


& Cordes (2003) quantify this effect by comparing the 
signihcance of detection through a FFT to that of single 
pulses, for any source with a given period and burst rate. 
They derive an analytic expression for the ratio of single¬ 
pulse to FFT detecti on significance, (S/N)sp/(S/ N)fft- 
Moreover, they plot ( McLaughlin fc Cordes|200^ Figure 
13) the regimes in which a single-pulse search is more sen¬ 
sitive than a FFT, which indeed shows that as a source 
emits more pulses, a FFT becomes more effective at de¬ 
tecting it. 

Similarly, [Keane] (|2010| uses this formalism to de¬ 
duce the regions of period/nulling fraction space in 
which single-pulse searches are more effective than FFT 
searches, and thus discoveries would likely be called 
RRATs. Using these constraints, he then finds that the 
average observed RRAT period is greater than the aver¬ 
age detected pulsar period. 

Due to this bias towards longer periods for RRATs, 
a statistical comparison of the periods of the detected 
RRAT and pulsar populations is not necessarily repre¬ 
sentative of the underlying distributions of these sources. 
However, we can characterize the detected populations 
and test whether those samples of periods are drawn from 
the same population. 

We use the two-sample Kolmogorov-Smirnov (K-S) 
test to compare the distributions of the properties men¬ 
tioned above. We add the sample of 20 RRATj^ found 
in this work to the 61 RRATs listed in the RRS^Talog. 
We then perform a K-S test between the periods of the 
RRATs and the ^2000 radio pulsars listed in the ATNF 
Catalogue (where we have excluded MSPs and magne- 
tars). This K-S test yields a p-value of 6 x 10“^^, imply¬ 
ing that the samples are drawn from extremely different 
distributions, as expected. 

We can also compare the periods of the new RRATs to 
those of all known RRATs, listed in the RRATalog. This 
yields a p-value of 0.16, meaning that the samples are 
consistent with coming from the same distribution. This 
is not surprising, since the sensitivities of recent surveys 
to RRATs are comparable to one a nother (e.g. Pulsar 


Arecibo L-band Feed Array survey, jPeneva et a. 


Parkes intermediate a nd high latitude surveys, 
Spolaor fc Bailes||2010 PMPS, Keane et al.||2Q10' 


Burke 


Finally, we compare the periods of the new RFtATs to 
those of the 128 pulsars discovered in the Drift-scan and 
GBNCC surveys. This yields a p-value of 0.002, imply¬ 
ing that the detected RRAT and pulsar period samples 
are likely drawn from different distributions. We plot 
the period distributions of the RRATs and pulsars found 
in each of the Drift-scan and GBNCC surveys in Fig¬ 
ure (left). These histograms have been normalized to 
account for the total number of sources. 

While we find that RRATs have preferentially longer 
periods than canonical pulsars (i.e. pulsars which were 
initially detected in a periodicity search), we mention 
above that this may be related to a classification bias. We 


We use 20 of the 21 RRATs presented here, since we exclude 
RRAT J0054+69 for which we could not reliably measure the pe¬ 
riod, as described in the text. 
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can avoid this bias by comparing the periods of single¬ 
pulse emitting pulsars to those of pulsars from which 
single pulses are not observed. Since single-pulse searches 
have uniform sensitivity regardless of the spin period of 
a pulsar/RRAT (as long as the pulses are identical), if 
fast rotators indeed emit bright RRAT-like pulses, we 
expect to detect them. We plot in Figure the period 
distribution of pulsars discovered in the GBNCC survey 
as well as known and newly discovered pulsars detected 
in the Drift-scan survey. We distinguish between sources 
from which single pulses were and were not detected. 
We see that generally, rapidly rotating pulsars do not 
emit bright single pulses, whereas slower pulsars often 
do. Moreover, a K-S test between these samples yields a 
p-value of 10“^^, indicating different distributions. 

Again, these distributions may be influenced by detec¬ 
tion biases. The S/N of a single-pulse detection depends 
on the width and amplitude of a given pulse. For a fixed 
flux density and duty cycle, we find that the flux per 
pulse is the same regardless of pulsar period, and there¬ 
fore expect only a mild bias in the distribution of sources 
detected via single pulses, due to broader pulse widths in 
slow pulses. On the other hand, if the d uty cycle varie s 
with period, for example as (e.g. Rankin 1993), 


then we expect longer-period pulsars to have larger flux 
per pulse, which would cause some bias in the observed 
distribution. Although these factors account for some of 
the difference between the period distributions, the ob¬ 
served effect (Figure is so pronounced that it could 
well also be caused by an intrinsic difference between 
the slow and fast pulsar populations and suggests that 
bright single-pulse emitters are preferentially slower ro¬ 
tators. Detailed quantitative follow-up with corrections 
for differences in pulse widths and fluences will allow firm 
conclusions. 


5.3. Dispersion measures 

We now turn to examining the biases affecting the DM 
distribution of a given survey’s discoveries. This prop¬ 
erty is heavily survey-dependent, since DM varies greatly 
depending on the surveyed region due to the non-uniform 
distribution of free electrons throughout the Galaxy. For 
this reason, we compare only the DMs of sources found 
in the same survey. 

Another effect that must be considered is the depen¬ 
dence of a survey’s sensitivity on DM. Since interstellar 
scattering timescales increase with DM an d with decreas - 


ing observing frequency (as / e.g. Rickett 199Q[ ) 


surveys are typically limited in their sensitivity to short- 
period pulsars at high DMs. This effect is mostly uniform 
for pulsars and RRATs, if their intrinsic period distri¬ 
butions are the same, although the narrower pulses of 
RRATs would likely be affected more by the diminished 
S/N due to scattering. 

The DM distributions of RRATs and pulsars discov¬ 
ered in each of the Drift-scan and GBNGG surveys are 
shown in Figure 12 (right). Note that these histograms 
have been normalized to account for the total number 
of sources. Performing a 2-sample K-S test between the 
DMs of the RRATs and pulsars in the Drift-scan survey, 
we obtain a p-value of 0.68, indicating that the samples 
are consistent with coming from the same DM distribu¬ 
tion. Repeating this for the sources in the GBNGG sur¬ 
vey, we obtain a p-value of 0.97, again implying a shared 


underlying DM distribution. 

5.4. Spatial distributions 

We can similarly compare the spatial distributions of 
RRATs and pulsars. Since different surveys cover differ¬ 
ent areas of the sky, meaningful comparisons are again 
only between sources found in the same survey. 

We compare the Galactic longitudes (/) and latitudes 
(h) of pulsars to those of RRATs, for the Drift-scan and 
GBNGG surveys. A K-S test for the I and b distri¬ 
butions in the Drift-scan survey yields p-values of 0.88 
for both, implying highly consistent spatial distributions. 
The same test for the GBNGG sources yields p-values of 
0.31 (/) and 0.40 (6), again implying that the two samples 
are consistent with having a shared spatial distribution. 

Examining the spatial distributions of sources is of in¬ 
terest since the Galactic position of a source may help 
constrain its age. For example, we expect young pulsars 
to be within the Galactic plane, where they formed, while 
older pulsars are found outside of the plane. The above 
results thus suggest that, according to spatial distribu¬ 
tions at least, RRATs and pulsars are consistent with 
having similar underlying properties. 

5.5. Fast radio bursts 


We did not find any FRBs in either the Drift-scan or 
the GBNGG survey. Since there are only several FRBs 


known (iThornton et al. 

2013 

Spitler et al. 

2014 Burke- 

bpolaor Bannister||20] 

L4|), their DM distribution is not 


well constrained. Moreover, tne publisned t rates 
have all been inferred from results of surveys conducted 
at a frequency of 1.4 GHz, making it difficult to predict 
a rate for surveys at 350 MHz since the spectral index 
for FRBs is not known. 

Assuming FRBs to have a uniform distribution of DMs, 
we can use the total observing time of each survey and 
telescope beam size to obtain an upper limit on the FRB 
rate, given that none have been detected. For the Drift- 
scan, assuming a Poisson process, we find a 99% upper 
limit on the FRB rate of ^ 1 x 10^ sky“^ day~^ above 
the minimum flux density, S rn\n = 260 mJy. For co m- 


parison, the rate inferred by Thornton et al. ( 2013) is 
1 X 10^ sky~^ day~^, where as the rate inferred by Burke- 


Spolaor & Bannister (2014) is 2 x 10^ sky“^ day“^, both 
for a huence of ^ 3 Jy s. For much of the GBNGG data 
analyzed with RRATtrap, the maximum DM searched was 
500 pccm“^, which is typic ally lower than the DM s of 
FRBs observed thus far (e.g. |Thornton et al.|2013|). We 
therefore do not quote here an event rate limit based on 
the GBNGG survey. A more thorough analysis will be 
presented in the future, placing better constraints on the 
FRB rates based on these surveys. 

6. CONCLUSIONS 


We have developed and used a new single-pulse sifting 
algorithm, RRATtrap, in order to find RRATs in GBT 
Drift-scan and GBNGG pulsar survey data, resulting 
in the discovery of 21 RRATs. We have characterized 
these RRATs by refining their periods, DMs, and burst 
rates. The new RRATs have DMs ranging from 15 to 
97 pccm“^, periods of 240 ms to 3.4 s, and burst rates 
of 20 to 400 pulses hr“^ at 350 MHz. We have also 
obtained phase-coherent timing solutions for four of the 
new RRATs. 
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Fig. 12. — Comparison of period (left) and DM (right) distributions of RRATs (blue-hatched) and pulsars (red) found in the Drift-scan 
and GBNCC surveys. The histograms hav e bee n n orm alized and include 11 Drift-scan and 10 GBNCC RRATs, and 35 Drift-scan and 93 
GBNCC pulsars. As described in Sections |5.2| and |5.3| the observed period distributions of pulsars and RRATs differ, whereas their DM 
distributions are similar. 



Fig. 13.— Period distribution for GBNCC discoveries and Drift- 
scan discoveries and detections of known pulsars. Sources de¬ 
tectable through single pulses are plotted in gray, with RRATs 
highlighted in dark gray. Sources not detectable through single 
pulses are outlined in black. All sources except the RRATs were 
detected in a periodicity search. 

We have presented our new RRAT discoveries in the 
context of the known RRAT population, as well as the 
radio pulsar population. We described the various detec¬ 
tion biases that exist between RRATs and pulsars, and 
took note of these when conducting statistical compar¬ 


isons of the two populations’ properties. We found that 
the DM and spatial distributions of RRATs and pulsars 
discovered in the same survey are consistent with having 
a shared underlying distribution. In contrast, the pe¬ 
riod distribution of detected RRATs is significantly dif¬ 
ferent than that of detected canonical pulsars. This may 
be an effect of the classification biases related to period 
and burst rate values of RRATs, and is not necessarily 
representative of the underlying RRAT and pulsar pop¬ 
ulations. We also find evidence that slower pulsars (i.e. 
P > 200 ms) are more likely to emit RRAT-like bright 
single pulses than are faster pulsars (P < 200 ms), al¬ 
though this conclusion is tentative. 

We have also shown that RRATs have a distribution of 
burst rates which can affect their classification as single¬ 
pulse or periodicity sources and also bears consequences 
for their discovery in surveys. We note that while we ob¬ 
served a certain distribution of burst rates, the underly¬ 
ing distribution is not known due to heavy observational 
biases against detecting RRATs with low burst rates. 
However, the increasing diversity in properties such as 
intermittency timescales of transients has begun to hint 
at a continuum of behavior, supporting the idea that 
RRATs are related to nulling pulsars and may simply 
represent extreme cases of the latter. 
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